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We predict that the maximal, steady-state ion capacity W(A) of radio-frequency (rf) traps, loaded 
at a rate of A particles per rf cycle, shows universal, nonlinear, nonmonotonic behavior as a function 
of loading rate A. The shape of Ns{X), characterized by four dynamical regimes, is universal, i.e., it 
is predicted to manifest itself in all types of rf traps independent of the details of their construction 
and independent of particle species loaded. For A <C 1 (Region I), as expected, Ns{\) increases 
monotonically with A. However, contrary to intuition, at intermediate A 1 (Region II), Ns{X) 
reaches a maximum, followed by a local minimum of Ns{X) (Region III). For A S> 1 (Region IV), 

Ns{X) again rises monotonically. In Region IV numerical simulations, analytical calculations, and 
experiments show Vs (A) ~ A^^®. We confirm our predictions both experimentally with MOT-loaded 
Na+ ions stored in a hybrid ion-neutral trap and numerically with the help of detailed ab-initio 
molecular-dynamics simulations. 

PACS numbers: 37.10.Ty, 52.27.Jt, 52.50.Qt 


I. INTRODUCTION 

Radio-frequency (rf) traps n are important devices 
in widespread use for the long-time storage of charged 
particles. These traps come in a multitude of shapes 
and sizes and their applications range from high- 

resolution spectroscopy Q to atomic clocks Q and quan¬ 
tum computers . They have also been used in nonlinear 
dynamics for the investigation of nonlinear phenomena 
ranging from crystallization [lol - [T5 | to the investigation 
of strange attractors [l^. A fundamental problem of 
great theoretical and practical interest is the maximal 
ion capacity Ns{X) of rf traps loaded at a constant rate 
of A particles per rf cycle. While the absolute value of 
Ns depends on the trap’s physical size and the details of 
its construction, the qualitative dependence of Ns{X) on 
A does not. In fact, we found that the shape of Ns{X) 
is universal, i.e., it is the same for any kind of rf trap, 
and shows four clearly defined dynamical regimes, which 
we label Regions I to IV. In this paper, based on physi¬ 
cal arguments and detailed ab-inito molecular-dynamics 
simulations, we predict the qualitative shape of the uni¬ 
versal curve Vs (A), and experimentally verify our pre¬ 
dictions with the help of a MOT-loaded ion-neutral trap 
[Tsil . While knowledge of the steady-state ion capacity Ng 
of rf traps in general is in itself an important fundamen¬ 
tal problem, the results are also of practical interest in 
atomic physics. Collision-rate experiments [Tsl - I^ . e.g., 
use the steady-state ion capacity to measure the total col¬ 
lision rate, because it ensures constant density, size, and 
temperature conditions during the measurement. This 
is particularly helpful when working with optically dark, 
closed-shell ions such as Rb+ and Na"*". 

Our paper is organized as follows. In Sec. [Ill we 
present the basic dynamical equations for the three- 
dimensional (3D) Paul trap and the linear Paul trap used 


in our molecular dynamics simulations together with the 
methodology according to which our simulations are per¬ 
formed. In this section we also present the theoretical ev¬ 
idence for the four nonlinear loading regimes encountered 
in these two trap types. In Sec. lIIIl we present our exper¬ 
imental evidence that confirms the prediction of the four 
different dynamical regimes. In particular, in the case of 
the linear Paul trap, we confirm the presence of the dip 
characterizing the dynamical Region III. In Sec. IIVI we 
discuss our results. We summarize and conclude our pa¬ 
per in Sec. s We also provide an appendix in which we 
derive the differential equation whose stationary solution 
yields the fundamental scaling relation Ns{X) ^ A^^^ in 
Region IV. 


II. THEORY 

Initially, we discovered the non-linear, nonmonotonic 
structure of Ns{X) in molecular-dynamics loading simu¬ 
lations of a 3D Paul trap. Denoting by rg and zq the dis¬ 
tances of the ring electrode and the end-cap electrodes 
from the trap’s center, respectively, by Uq and Vg the 
amplitudes of the dc and ac voltages applied to the trap, 
respectively, by P the damping constant (generated, e.g., 
by laser cooling M), and by oj the (angular) frequency 
of the trap’s ac voltage, the (dimensionless) equations 
of motion of N particles in the trap [r = {x,y,z)] are 
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where i = 1,..., N labels the trapped particles, 
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are the two control parameters of the Paul trap 
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where Q is the charge and m is the mass of each of the 
trapped particles, cq is the permittivity of the vacuum, 
and 

2r 

7 = roP = — (5) 
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is the dimensionless damping constant. We use dimen¬ 
sionless quantities in this section because only this way 
is it possible to see that the equations of motion of par¬ 
ticles in the Paul trap do not depend on the six phys¬ 
ical parameters rg, zq, P, oj, Q, and m separately, but 
only on the three scaled, dimensionless parameters a, q, 
and 7 - Therefore, instead of the need to explore a six¬ 
dimensional parameter space, which is practically impos¬ 
sible, we only need to explore a three-dimensional param¬ 
eter space. This grows to four dimensions, if we include 
the loading rate A. 

In our simulations particles are created, one at a time, 
at times k = 1 , 2 ,..., either with zero initial velocity 
(a good approximation for MOT-loaded ions [1^) or with 
a thermal velocity distribution as discussed below. As¬ 
suming that the creation times are uncorrelated, the time 
intervals Atk = —tk are Poissonian distributed with 
probability distribution P{At) = Aexp(—AAt). Con¬ 
cerning their spatial distribution, we assume that the 
particles are created at random positions with uniform 
distribution within an ellipsoidal volume with semi-major 
axes Lx, Ly, Lz, centered at the origin of the trap. This 
includes the case of a spherical loading zone of radius 
i?, in which case we have R = Lx = Ly = Lz- In the 
time interval between any two creation events, i.e., for 
tk < t < tk+i, the particles in the trap are governed by 
the equations of motion O- When arriving at tfe+i, and 
before creating the next particle, we check whether one 
or more particles have left the trap by crossing an ab¬ 
sorbing boundary. The existence of an absorbing bound¬ 
ary is a fundamental property of all rf traps, which de¬ 
termines and limits the storage capacity of any given 
trap. This boundary may be due to any number of un¬ 
avoidable physical causes, such as the trap’s electrodes, 
or instabilities induced by higher-order rf multipoles 
Since all traps are constructed differently, and to show 
that our predicted effect is robust with respect to vari¬ 
ous geometries of absorbing boundaries, we used spher¬ 
ical boundaries of radius i?sph, boxes with side lengths 


FIG. 1: (Color online) Number of particles N{t) in a 3D Paul 
trap as a function of t (in rf cycles), loaded with loading rate 
A = 1 particle / rf cycle inside of a spherical volume of ra¬ 
dius i? = 3 centered at the origin of the trap. Particles are 
absorbed at a critical distance Rc = 15. Trap control pa¬ 
rameters are 'y — 0, q = 0.2, a = 0.02. Inset: Test of the 
loading equation for a loading rate A = 300 particles / rf 
cycle. Red solid line: Number of particles N(t) obtained via 
numerical simulation of the loading process. Blue dotted line: 
Prediction according to ®. 


2xboxj 2yboxj 22boxj and cylinders with radius i?cyi in the 
x-y plane and length 22cyi = {[q^/{^B)] — 1/2}^/^ in 2 
direction to cover a wide variety of possible boundary 
geometries. Following instantaneous deletion of all parti¬ 
cles that exceed the confines of the absorbing boundary, 
the next particle is loaded at f = tk+i- This procedure 
is followed for all tk until a pre-specified maximal simu¬ 
lation time is reached. 

The result of a typical trap loading simulation for a 
spherical loading zone with i? = 3 is shown in Fig.[T] The 
red, fluctuating line in Fig. [T] shows the time evolution 
of the particle number N{t) in the trap for A = 1 as a 
function of t (in rf cycles) for q = 0.2, a = 0.02, and 
.Rsph = 15. A near linear rise of N{t) is followed by a 
sharp bend into a steady state in which N{t) fluctuates 
around Ng = {N{t))t, the time average of N{t) in the 
steady-state. 

For many applications (see, e.g., [11,0) it is necessary 
to know Ns{X) for a given rf trap. Figure [2] shows W(A) 
as a function of loading rate A for the same Paul trap 
control parameters used to generate Fig. [T] The result¬ 
ing curve in Fig. [5] clearly shows four distinct regimes, 
labeled I to IV. Region IV is the most straightforward 
to understand physically; therefore, we discuss it first. 
In Region IV, the loading rate is so large that a large 
space-charge density develops in the loading region. The 
resulting large electric field accelerates the loaded parti¬ 
cles outward toward the absorbing boundary where the 
particles are lost from the trap. Since in this case the 
forces due to the space charge completely overwhelm the 
forces due to the trap fields, the particles’ d yna mics are 
accurately described as a Coulomb explosion . Denot¬ 
ing the number of particles inside of the loading zone by 
N{t) [in contrast to N{t), which refers to the total num- 
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FIG. 2: (Color online) Nonmonotonic loading curve Ns{X) for 
a 3D Paul trap with 'y — 0, q — 0.2, a = 0.02, R = S, and 
I?sph = 15. Solid green dots: Results of 3D molecular dynam¬ 
ics loading simulations. The red bars indicate the amplitudes 
of the N{t) fluctuations in the saturated state. For the first 
two dots the N{t) fluctuations are smaller than the plot sym¬ 
bols. The four distinct regions of the loading curve are labeled 
I to IV. The red solid line is the curve Ns{\) = 5.2 x 
which confirms the ~ behavior of 3D Paul traps in Re¬ 
gion IV. 



FIG. 3: (Color online) 3D Paul trap capacities Ns{X) for var¬ 
ious trap control parameters q and 7 = 0. Solid red squares: 
q — 0.1; open blue circles: q = 0.15; solid cyan circles: 
q = 0.2; open black triangles: q = 0.25; solid grey trian¬ 
gles: q = 0.3. The trap parameter a associated with each of 
the five q values is a — q^/2. The four dynamical regimes 
are labeled I to IV. The solid red lines connecting the data 
points are drawn to guide the eye. All five curves asymptote 
to Nb{X) ~ A^^® for large A. 


ber of particles in the trap], the temporal evolution of 
the number of particles, N(t), inside of the loading zone 
is governed by the differential equation (see appendix): 




( 6 ) 


where R = [16/(97r)]^/^i? Ri R is the effective radius 
of the loading volume and A is the number of parti¬ 
cles loaded per unit of dimensionless time. According 
to (|A.1I) . A is related to A, the number of particles loaded 


per rf cycle by A = A/tt. In the stationary state, we have 
dN{t)/dt = 0. Therefore, we obtain from (jS]) 

iV,(A) = i?A2/3= j i?A2/3. (7) 

Because of continuity, the steady-state number of par¬ 
ticles in the trap, Ns{X), is proportional to the number 
of particles in the loading region, Ns{X). Therefore, the 
A^/^ dependence of Ns{X) is reflected in Fig. [5] (solid red 
line). Because there are more particles in the trap than 
there are in the loading zone, the pre-factor 5.2 of A^/^, 
stated in the caption of Fig. [21 is larger than the pre¬ 
factor r: 1 of A^/^, computed from (O with 

R = 3. 

The solution of ([6]) can be stated implicitly in closed 
form: 
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where a = That (jS)) is indeed a solution of (jH) 

may be checked immediately by differentiating ([5]) with 
respect to t. Since the trap potentials are not important 
in Region IV, the results ([6]) - dH) apply universally to 
all rf traps, for instance 3D or linear Paul traps. 

In order to test (jS]) , we chose R — 3 and a large loading 
rate of A = 300 particles / rf cycle to be sure that we are 
in Region IV. The result of the corresponding molecular 
dynamics simulation of N{t) is shown as the solid, fluctu¬ 
ating red line in the inset of Fig. |TJ The blue, dotted line 
in the inset of Fig. |T] is the prediction according to ([5]). 
Both agree perfectly within the expected fluctuations of 
N{t), which are due to the Poissonian loading process. 

While for large A the trap potentials are not impor¬ 
tant, they become progressively more important when 
A is lowered. In this case, for low enough A, particles 
created close to the edge of the loading zone no longer 
have enough energy to overcome the trap potentials and 
are reflected back into the interior of the trap. When 
back-reflection occurs, fewer particles escape and Ns{X) 
increases for decreasing A up to a maximum in Region 
II, effectively creating a dip in Region III (see Fig. |2]). 
However, the reflected particles will not stay in the trap 
forever. Due to rf heating [l^, [2l|, and given enough 
time (small loading rates), these particles will heat out 
of the trap, eventually lowering the number Ns{X) of 
stationary-state particles in the trap below the bottom 
of the dip in Region III (see Fig.|2]), thus explaining both 
the formation of the maximum in Ns{X) (Region II in 
Fig. 12]) and the eventual decline of Ns{X) in the direction 
of ever diminishing loading rates (Region I in Fig. |2]). 
Since rf heating is a universal feature of all rf charged- 
particle traps, we predict that the qualitative shape of 
Ns{X) is universal for all rf traps. 
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FIG. 4: (Color online) Saturation curves Ns{X), for a linear 
Paul trap. Red pluses: q — 0.35; blue crosses q = 0.40; cyan 
diamonds: q = 0.45; black open squares: q = 0.50; grey solid 
squares: q — 0.55; green open circles: q = 0.60; magenta solid 
circles: q — 0.65; brown open triangles: q — 0.70. The solid 
red lines connecting the data points are drawn to guide the 
eye. All eight curves asymptote to Ns{X) ~ for large 
A. The labels I to IV refer to the four different dynamical 
regimes. 


To strengthen the claim of universality of the non¬ 
monotonic curve shown in Fig. [5J we show in Fig. [3] 
the saturated number of particles Ns{X) as a function 
of loading rate A for several values of q with a = q'^/2, 
Lx — Ly = 3, Lz — 6, and Xuox — ^box — -^^box — 15. 
All five curves in Fig. [3] clearly show all four dynami¬ 
cal regimes, again labeled I to IV. This shows that the 
nonmonotonic behavior is robust with respect to (a) a 
change in q, (b) a change in the geometry of the load¬ 
ing zone (spherical in Fig. [31 elliptical in Fig. [3]), and 
(c) a change in the geometry of the absorbing boundary 
(spherical in Fig. |2l cubic in Fig. [3]). This applies in par¬ 
ticular to Region IV, in which all five curves in Fig. [3] 
are seen to converge to the same ^ X^/^ asymptote. The 
independence of a and q in region IV is explained by ([HI) , 
which is independent of the trap potentials and therefore 
independent of the trap parameters a and q. Figure [3] 
also shows that the dip becomes shallower with both in¬ 
creasing and decreasing q and shifts to the right with 
increasing q. In addition, we see that Region II becomes 
more extended with decreasing q. At present we do not 
have a theoretical explanation for these observed effects. 

To check the universality of the nonmonotonic behav¬ 
ior, we also simulated a linear Paul trap [2^ whose equa¬ 
tions of motion, in the notation and units of o, are given 

by 


Xi + la — 2flsin(2t)la;i — BxA ^ ^ 

iji-[a- 2qsm{2t)]y, - By A = ^ (9) 

Zi-\-2Bzi j j=i T* 


where 5 is a positive constant. We simulate the loading 
process of the linear Paul trap in analogy to the 3D Paul 
trap as discussed above with a loading region of radius 
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FIG. 5: (Golor online) Experimental saturation curve for the 
hybrid trap El. operated with q = 0.26. All four predicted 
dynamical regimes are present. The red solid line is the curve 
Ns = 48000A^^®. The error bars represent the statistical er¬ 
rors associated with multiple data runs. The zero-point of the 
horizontal scale is suppressed. 



i? = 3, located at the geometric center of the linear trap 
and a cylindrical absorbing boundary with i?cyi = 15. 
This time, however, we use the equations of motion (jH]) 
between creation times. Figure |3| shows the resulting 
Ns{X) for eight different q values with o = 0, and B = 
0.042. We clearly see the four different regions of Ns(X), 
observed previously in the case of the 3D Paul trap. This 
provides corroborating evidence for the universal nature 
of the shape of Ns{X) for all types of rf traps. 


III. EXPERIMENT 


The ultimate test of our theoretical predictions is an 
experiment. For this purpose we used a MOT-loaded 
ion-neutral hybrid trap ES- Figure |S] shows the results 
of our experiments. All four predicted dynamical regimes 
are present and the shape of the experimental Ns{X) 
curve is seen to be qualitatively the same as predicted 
by the model simulations and our qualitative analysis of 
the physical mechanisms that determine the steady-state 
populations in rf traps. From the experimental data we 
obtain e = 0.745±0.098 for the exponent e in Ns{X) ^ A*^, 
which is consistent with the predicted value e = 2/3 in 
Region IV (red solid line in Fig. |S]). We note that, cor¬ 
roborating the universality claim, the qualitative shape 
of the nonlinear, nonmonotonic behavior of Ns{X) was 
observed in all of our simulations and experiments, inde¬ 
pendent of rf amplitudes, rf frequencies, and MOT sizes. 
With a modulation depth of about 5%, the dip in Fig. |5| 
is small. However, noticing that the experimental q of 
q = 0.26 is small, this observation is consistent with 
Fig. m which shows that the modulation depth of the 
dip decreases with decreasing q. 
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IV. DISCUSSION 

To strengthen the universality claim and to emphasize 
the robustness of our predictions, we performed numer¬ 
ous additional simulations that all confirmed that the 
qualitative shape of the saturation curve, characterized 
by its four dynamical regimes, is insensitive to both the 
type of traps used and their particular loading mecha¬ 
nisms. In particular, we performed additional simula¬ 
tions with (a) various shapes of the loading ellipse with 
aspect ratios of up to 1 : 1 : 4, (b) re¬ 

placing the Poissonian distribution of loading times with 
a uniform distribution, (c) replacing the uniform spatial 
distribution with a Gaussian distribution, (d) various ge¬ 
ometries of the absorbing boundary and boundary loca¬ 
tions, and (e) increasing the number of trapped particles 
to up to 1000 by changing the diameter of the absorbing 
boundary. All five numerical tests confirmed the qualita¬ 
tive shape of Ns{X) as shown in Figs.[2]-[5] In addition, 
we performed the following two checks concerning ion 
cooling and the effect of temperature. 

In many ion-trap experiments, strong laser cooling is 
continuously switched on, even during the loading stage 
(see, e.g., [l2|)- Therefore, to test the influence of cool¬ 
ing during trap loading, we ran additional 3D Paul-trap 
loading simulations with strong damping switched on, 
so strong in fact, that it would crystallize the ions [ 13 - 
[ 12 , [13 if we were not constantly loading. This corre¬ 
sponds to 7 values between 5 x 10“^ and 10“^. As a re¬ 
sult we find that we still have all four dynamical regimes, 
in particular the dip. Moreover, the Ns{X) curves in the 
two cases (with and without damping, respectively) are 
nearly identical, differing from each other only within the 
natural fluctuations of Ng due to the Poissonian loading 
process. In some of our more recent simulations, just 
to push the envelope, we increased the damping to five 
times the one needed for crystallization (7 = 5 x 10“^). 
It still did not have any effect on the qualitative shape of 
the loading curves. Thus, we conclude that the presence 
of laser cooling that can realistically be achieved experi¬ 
mentally has no effect on our predicted phenomena. We 
also mention that a large laser cooling power does not 
necessarily result in a large 7 , since 7 represents the bal¬ 
ance between laser cooling and the substantial amount 
of heating caused by the loading process. 

In our experiments we load from a cold MOT and 
the assumption of zero kinetic energy at the instance of 
charged-particle creation is justified. However, in case 
the trap is loaded by ionizing the rest gas or via a ther¬ 
mal neutral beam, the particles may be created with sub¬ 
stantial initial kinetic energy. To test the effect of initial 
kinetic energy on our predictions, we performed addi¬ 
tional loading simulations, imparting a random velocity 
on each particle at the instance of its creation whose en¬ 
ergy equivalent was up to 20% of the trap depth. In 
our experiments this would be equivalent to about room 
temperature. Even in this case we are still able to ob¬ 
serve all four dynamical regions. The explanation is that 


in steady-state there is an equilibrium between particles 
created and particles leaving the trap. Those that leave, 
already have kinetic energies of the order of the trap’s 
depth, which they are able to impart, via Coulomb colli¬ 
sions, to newly created particles. Therefore, the fact that 
particles are created with a kinetic energy less than the 
trap’s depth is only a minor perturbation on the ener¬ 
getic particle dynamics that is already going on inside of 
the trap. 

In our simulations we found that the dip in Region III 
is more pronounced if the absorbing boundaries are fur¬ 
ther from the edges of the loading zone. This observa¬ 
tion explains why the dip is more pronounced in Fig. [2] 
compared with the dips in Fig. [3] In Fig. |3] we used 
an elliptical loading zone with = 6 compared with 
the spherical loading zone in Fig. [2] with i? = 3. There¬ 
fore, in z direction, it is easier for the particles in Fig. [3] 
to bridge the gap to the absorbing boundary, because 
many of them are already created closer to the absorb¬ 
ing boundary. The shorter distance to cover results in a 
shallower dip. 


V. SUMMARY AND CONCLUSIONS 

In this paper we report the discovery of the non¬ 
monotonic shape of the saturated ion number iVs(A) of 
rf traps as a function of loading rate A and present ev¬ 
idence for its universality. Four dynamical regions are 
predicted. In Region I, Ns{X) increases monotonically 
with A, reaching a maximum (Region II) at intermediate 
loading rates A, followed by a valley (Region III), and an 
ultimate ^ A^/^ increase of Ns{X) for very large loading 
rates. We argue that the four regions are expected on 
the basis of physical reasons and are caused by the inter¬ 
play between trap potentials, space-charge effects, and 
rf heating. The validity of our predictions, in particular 
their universality, are corroborated with the help of ab- 
initio molecular-dynamics simulations of a 3D Paul trap 
and a linear Paul trap, which both show the predicted 
qualitative shape of Ns{X). The theoretical predictions 
are confirmed experimentally with the help of Na'*" ions 
in a MOT-loaded linear Paul trap. 

Apart from our research group , several other 
groups [l 6 l - [^ have the necessary experimental facilities 
to test our predictions in the case of MOT-loaded traps, 
and since, as we showed, the predicted phenomenon is ro¬ 
bust with respect to rf trap types, loading mechanisms, 
and temperature effects, we hope that other research 
groups may soon test and confirm our predictions. 
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Appendix: Region IV differential equation 


The general solution of (IA.5I) is 


In this appendix we present a simple, explicitly solv¬ 
able, analytical model that reproduces the scaling 
of Ng in Region IV, the fast-loading regime. While in 
Sec.ini we used dimensionless quantities to bring out the 
scaling properties of the Paul-trap equations, it is more 
convenient in this appendix to derive our equations us¬ 
ing SI units. In these units, we denote time (measured 
in seconds) by t, the loading rate (measured in parti¬ 
cles per second) by A, and the radius of the loading zone 
(measured in meters) by R. However, in order to make 
contact with the formulas in Sec. m it is convenient to 
define A, as we did in Sec. m as the number of particles 
loaded per rf cycle, and by 

A = tqA = 2A/tLi = A/tt (A.l) 

the dimensionless loading rate per unit of dimensionless 
time. The unit of time, tq, in (lA.ll) is defined in (I3|). The 
last equality in (lA.ll) comes about since in dimensionless 
time an rf cycle has a length of tt [see Cl)]- 

We are now ready to start our derivation of equation 
® in Sec. IE which holds in Regime IV. As discussed 
in Sec. IE in this regime we may neglect the trap poten¬ 
tial altogether. We assume that the trapped particles are 
created at random locations inside of a sphere of radius 
R with uniform probability distribution and loading rate 
A. If there are N particles present inside the sphere of ra¬ 
dius R, the charge density, approximated as a continuous 
distribution, is 


r(T) = Aexp(\/i^r2r)-|-Rexp(—V^rir), (A.7) 

where A and B are constants. If we assume that at r = 0 
the particle is created a distance s away from the center 
of the loading sphere, (IA.7|) may be written as 

r(T) = s cos h(V^Hr), (A.8) 

and the velocity of the particle is 

f(r) = sinh(\/i^r2T). (A.9) 

Define Tq as the time it takes the particle to reach the 
rim of the loading sphere at radius R if the ion starts at 


radius s with zero velocity. Then; 

R = scosh{\/^QTo) (A.10) 

^ cosWV^DTo) = -. (A.ll) 

s 

When it arrives at i?, the velocity of the ion is 

Vo = s^/nQ sinh( V^DTo). (A.12) 

Use (lA.llll . together with cosh^(a:) — sinh^(a:) = 1, to 
write (IA.12I) in the form 

Vo = ( a . 13 ) 


3iVQ 

4:TTR^ 


(A.2) 


where Q is the charge of each trapped particle. Using 
Gauss’ law, the radial electric field, pointing outward, a 
distance r away from the center of the sphere, is 


Next, we compute the average velocity vo with which a 
randomly created particle arrives at R. Denoting by V 
the volume of the loading sphere, V = AirlR/S, and using 
the fact that we assume a uniform probability distribu¬ 
tion of particle creation positions within the sphere, we 
obtain: 


E = 


pr 

Sen 


(A.3) 


Therefore, the radial, outward directed force experienced 
by an ion a distance r away from the loading sphere is 


F = 



( NQ^ \ 

y47reo.R^/ 


(A.4) 


This equation shows that the force experienced by a sin¬ 
gle trapped particle is like the harmonic force of an in¬ 
verted oscillator. Therefore the equation of motion of the 
ion is 


mr = F 


f = Nfl^r, 


(A.5) 


% = / \/r^-s^ dV 

^ Jv 

R3 Jo ^ 16 

(A.14) 

Since the average ion arrives at R with an average ve¬ 
locity Vo, directed radially outward, the average number 
dNi of particles lost from the loading sphere in time dr is 
the number of particles in a shell of radius R and width 
vodr. Explicitly; 

N^/^ndr. (A.15) 


dNi = -^AnR^vodT = 

(qj 


Ott 

16 


where m is its mass and 


n = 



1/2 


The number of particles gained due to loading with rate 
A is 


dNg = Adr. 


47reomI?3 


(A.6) 


(A.16) 
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Therefore the total change dN in the number of particles 
present in the loading sphere is 

dN = dNg-dNi= KdT-CN^/^dT, (A.17) 

where 

C = (A.18) 

16 VdTreoTOi?^/ 

Dividing (IA.17I) by dr, we obtain a first-order differential 
equation for the number of particles N inside the loading 
sphere: 


dN , 

— = k-CN^/^. 

dr 


(A.19) 


To transform this equation into its dimensionless form, 
we use the unit of time, tq, and the unit of length, Iq, 


defined in m and o, respectively, which relates t. A, 
and R in SI units to their dimensionless counterparts, t, 
A, and R, respectively, according to [see also (lA.lIl 


T = tTQ, A = X/to, R=RIq. 
Using (IA.20I) in (IA.19I) and defining 


R = 



we arrive at 


dN (N{t)V'^ 


which is identical with (jb]) of Sec. m 


(A.20) 


(A.21) 


(A.22) 
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